Critical load and congestion instabilities in scale-free networks 
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We study the tolerance to congestion failures in communication networks with scale-free topology. 
The traffic load carried by each damaged element in the network must be partly or totally redis- 
tributed among the remaining elements. Overloaded elements might fail on their turn, triggering the 
occurrence of failure cascades able to isolate large parts of the network. We find a critical traffic load 
above which the probability of massive traffic congestions destroying the network communication 
capabilities is finite. 
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Complex heterogeneous connectivity patterns have 
been recently identified in several natural and technologi- 
cal networks Jl], ||, || . The Internet and the World-Wide- 
Web (WWW) networks, where nodes represent routers 
or web pages and edges physical connections or hyper- 
links, appear to have a topology characterized by the 
presence of "hubs" with many connections to peripherical 
nodes. Empirical evidence recently collected shows that 
this distinctive feature finds its statistical characteriza- 
tion in the presence of heavy-tailed degree distributions 

i> i & i H 0- In thc Intcrnet > for instance, 
the statistical analysis reveals that the degree distribu- 
tion .P(fc), defined as the probability that any node has 
k links to other nodes, is well approximated by a power- 
law behavior P(k) ~ fc"T, with 7 « 2.2 @, 0, This 
makes the Internet a capital example of the recently iden- 
tified class of scale-free (SF) networks Q. The statistical 
physics approach has been proved to be a valuable tool 
for the study of complex networks, and several interesting 
results concerning dynamical processes taking place on 
complex networks have been recently reported. In partic- 
ular, the absence of the percolation Jl2|, |l3) and epidemic 
fli] , [[5) thresholds in SF networks has a large impact 
because of its potential practical implications. The ab- 
sence of the percolating threshold, indeed, prompts to an 
exceptional tolerance to random damages fl6| . This is a 
property that assumes a great importance in communica- 
tion networks, guaranteeing the connectivity capabilities 
of the system. 

Percolation properties of SF networks refer only to the 
static topological connectivity properties (lj, |l3| . On the 
other hand, in the Internet and other communications 
networks, many instabilities are due to traffic load con- 
gestions |l7], The traffic load carried on the failing 
nodes or connections is automatically diverted to alter- 
native paths on the networks and instabilities can spread 
from node to node by an avalanche of traffic congestions 
and overloads. For instance, route flaps have led to the 
transient loss of connectivity for large portions of the In- 



ternet. These instabilities are thus of a dynamical nature 
and depend on how information is routed and distributed 
in the network. The models proposed so far, however, 
deal with regular structures |n| 20 ; EH) and do not take 



into account the complex topology of SF networks. 

In this paper, we propose a simple model aimed at the 
study of failure cascades generated by the redistribution 
of traffic load by congested links or nodes in SF networks. 
We find that the system behavior depends on the aver- 
age traffic load imposed to the network. Above a critical 
value of the average traffic load, a single failure has a 
finite probability of triggering a congestion avalanche af- 
fecting a macroscopic part of the network. The present 
analysis thus reveals the existence of a transition from a 
free phase to a congested one as a function of the amount 
of traffic carried by the network. Contrary to what hap- 
pens for the static percolation transition |l2|, loaded 
SF networks exhibit a finite threshold above which the 
system can develop macroscopic instabilities with respect 
to small damages if we consider the dynamics of the traf- 
fic carried on top of them. The results provided here 
represent a first step towards a more complete modeling 
of traffic instabilities in real communication networks. 

In order to include the degree fluctuations of SF net- 
works we shall use in the following the Barabasi- Albert 
model p2fl . This is a stochastic growth model in which 
one starts from a small number m of nodes and at each 
time step a new node is introduced. The new node is con- 
nected preferentially to m old ones (for the simulations 
we used m = 3) with a probability Il(fej) = ki/J2j kj pro- 
portional to the degrees ki of the nodes. The repeated it- 
eration of this scheme gives as a result a complex network 
with a topological structure characterized by a power- 
law degree distribution P(k) — 2m 2 fc~ 7 with 7 = 3 and 
average degree (k) = 2m. In principle, one might also 
consider a more general class of complex networks with 
variable power-law degree distributions ^ |). 

To simulate the flow of data packets on SF networks, 
taking into account the load redistribution in case of 
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damages, we need to specify the initial state of the net- 
work; i.e. the load of traffic flowing through each link. 
An estimate of such load, assuming that the routing takes 
place following the minimum path, is given by the total 
number of shortest paths between any two nodes in the 
network that pass through the node i. This magnitude is 
called betweenness or load [^3, 24 and has been recently 
studied in SF networks. This property of the network 
can also be defined in terms of links. In real systems, 
however, the amount of traffic carried by each link is 
a fluctuating quantity that depends on many variables 
such as number of users, routing agreements, and avail- 
able bandwidth. For this reason, we associate to each link 
connecting the nodes i and j of the network a load l^j 
drawn from a probability distribution that specifies the 
initial traffic load of the system. For simplicity, we have 
considered a uniform distribution U(£) for < £ < 1, 
taking the form 



U(£) 



^y, £e[0,2(e)} if <f) < 0.5 
jp^j, te [2(g) -1,1] if W > 0.5 



This uniform distribution implies that the minimum ini- 
tial load carried by a link is bounded by a nonzero value 
for an average load (I) > 0.5, which means that there 
will be no links with load smaller than this lower bound. 
The results reported in this paper were obtained with the 
initial distribution U{£). In order to test the universal- 
ity of the critical behavior, we have also considered the 
distribution F{£) = ({t)- 1 - 1)(1 - ^CW" 1 - 2 ), £ e [0, 1], 
which allows the existence of links with a very small load, 
irrespective of the average load (£) flowing through the 
system. Both distributions U(£) and F(£) yield the same 
qualitative results. Along with the load, we associate to 
each link the same capacity C that, without loss of gener- 
ality, we fix equal to one. This choice can be considered 
as a first approximation since the actual difference be- 
tween line bandwidths in communication networks can 
be large. In this perspective, we consider as the most 
important source of heterogeneity the flow of different 
amounts of load through the network. 

The dynamics of the model is defined by a simple 
threshold process. A link is selected at random and over- 
loaded by raising its traffic. When the load carried by 
a link is t^j > C, i.e. when it exceeds the link's ca- 
pacity, the link is considered congested and the load it 
carries is diverted among its (not overloaded) neighbor- 
ing links. This amounts to consider that the time scale of 
the local congestion is greater than the time scale charac- 
terizing the reorganization of the routing procedure. The 
redistribution of the load on its turn might provoke that 
other links become overloaded, thus triggering a cascade 
of failures. We have explored two physically different set- 
tings of the load redistribution rule. The first consists on 
equally distributing the load of a congested link among 




FIG. 1: Phase diagram of the system. Pg measures the prob- 
ability of having a giant component of active nodes Q of the 
order of the system size. Three different definitions of the 
model are represented: (a) random and dissipative, (b) ran- 
dom and conserved, and (c) deterministic and conserved. In 
all cases the size of the system has been set to N = 10 4 nodes 
(3 x 10 4 links). 



the non-congested neighboring links. We refer to it as 
the deterministic redistribution rule, respectively. The 
second case will be called random redistribution because 
when a link is overloaded, a random amount of load is 
redistributed to each of the remaining working links in its 
neighborhood. Finally we note that in the rare event in 
which the congested link has no active neighbors, its load 
can be equally shared among all the remaining working 
lines of the network or just be considered as lost from the 
network. This amounts to a conserved or dissipative re- 
distribution rule. Many physical systems display critical- 
ity only when energy is conserved ]25|, ^6) . In distributed 
networks such as the Internet, however, it is common to 
discard packets if there is not a route available at the 
moment. As we shall see in the following, the results do 
not depend qualitatively on the conserved nature of the 
traffic load. 

We have performed large-scale numerical simulations 
by applying repeatedly the rules stated above on BA 
networks. The sizes of the networks used in the simu- 
lations range from N = 5 x 10 3 nodes (15 x 10 3 links) 
to N = 10 5 nodes (3 x 10 5 links). All numerical re- 
sults have been obtained by averaging over 10 different 
networks and, at least, 100 different realizations of the 
initial load distribution. 

In order to inspect the occurrence of dynamic instabil- 
ities, we construct the phase diagram of the system. The 
order parameter can be identified as the probability Pg 
of having a giant component Q of connected nodes with 
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size of the order of the network size. The giant compo- 
nent is defined as the largest component of the network 
made by nodes connected by active links, after the sys- 
tem has reached a stable state (when I^a < C for all i 
and j). The existence of a giant component implies that 
a macroscopic part of the network is still functional. If 
the giant component of the network is zero, the commu- 
nication capabilities of the network are destroyed and a 
congestion of the order of the system size builds up. It 
is worth noticing that although we have determined the 
giant component size in terms of nodes and the dynami- 
cal rules of the model are expressed in terms of links, the 
results are completely equivalent since a connected node 
is defined as a node with at least one active link. 

In Fig. [j] we plot the order parameter Pg as a function 
of the average load (£) . At low values of the average load, 
the network always reaches a stable state in which the 
number of isolated nodes is very small and with prob- 
ability Pg — 1 the network has a giant component of 
connected nodes of the order of the system size. When 
increasing the load imposed on the network, the system 
starts to develop instabilities. In particular, above a crit- 
ical load (£)c, whose value depends on the model consid- 
ered, with a finite probability the system evolves to a 
congested state without giant component of connected 
nodes; i.e. the largest set of connected active nodes has 
a density of order TV -1 . This implies a probability of 
having a giant component Pg < 1, which is decreasing as 
the load is progressively increased. At an average load 
(I)*. 1 ~ 0.82 we get that Pg = 0, signalling that, with 
probability one, any instability will propagate until the 
complete fragmentation of the network. It is worth re- 
marking that this scenario is rather different from the 
percolation one in which the probability of having a gi- 
ant component is abruptly dropping from one to zero at 
the transition point. Here, the probability decays con- 
tinuously to zero and we have a wide region of {£) where 
the initial instability can trigger a destructive congestion 
with probability 1 — Pg. Fig. || illustrates the probabil- 
ity p{S) that the isolated network has a size S in the 
case that no giant component of connected nodes has 
survived. The distribution is rather peaked also at rela- 
tively small values of (£), almost affecting the totality of 
the network. 

The phase diagram obtained in Fig. [I] points out that 
the value of the average load at which Pg = is rela- 
tively high. On the other hand, the value at which Pg 
is appreciably smaller than one is well below the theo- 
retical capacity of the network measured as the capacity 
C = 1 of the individual links {(£)[ ~ 0.15 - 0.3, see 
Fig. |]) . This evidence defines a wide region of load val- 
ues in which there is a small but finite probability that 
a small instability propagates through the system and 
could explain why in real communication networks con- 
gestions of diverse degree can be encountered from time 
to time. In Fig. R] we also report curves obtained for the 
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FIG. 2: Distribution of the isolated network size p(S) as a 
function of S for several values of {£). The network is formed 
by N = 10 4 nodes. 
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FIG. 3: Cumulative distributions of avalanche sizes for dif- 
ferent values of the average load handled by the network 
for the random and dissipative definition of the model. The 
straight line in a linear-log plot indicates that the probability 
avalanche distribution follows a power law with exponent — 1 . 
The inset shows the scaling of the cumulative size of congested 
lines as the network grows in size for (I) = 0.25. 

different dynamical rules defined previously and it is in- 
teresting to remark that the highest level of stability has 
been obtained for the deterministic and conservative ver- 
sion of the model — 0.3). Moreover, the dynamics 
of the transition does not depend on the distribution of 
initial loads and the network size N . 

Another way to shed light on the congestion dynam- 
ics is to inspect the process of generation of the con- 
gestion instability. We define the size s of a congestion 
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burst or avalanche as the total number of simultaneously 
overloaded links. The cumulative distribution P(s) of 
avalanches of size larger than s for several values of the 
average load imposed on the network and four different 
system sizes have been plotted in Fig. || for the random 
and dissipative definition of the present model. The main 
plot is in log-linear scale, so that a straight line corre- 
sponds to a power law of the form p(s) ~ s _1 for the 
probability of observing an avalanche of size s. Power- 
laws with exponent — 1 have been found for several char- 
acteristic features of Internet traffic such as latency times, 
queue lengths, and congestion lengths ^7], [2£|. In 
the figure we focus on region close to the stable region 
(£) ~ 0.20, which means that the power-law behavior ex- 
tends to values far from the instability transition. This 
fact confirms that it is not necessary that the network 
operates very close to a critical point in order to observe 
power-laws in the distribution of several quantities. The 
inset in Fig. || shows that the cumulative size of over- 
loaded links also scales with the system size, the scaling 
dynamics, however, remaining the same. This may help 
understand why power-law distributions observed in real 
communication networks have been measured for differ- 
ent network sizes, i.e. both for local networks and for 
networks that extend to a very large scale. 

In summary, we have introduced a simple threshold 
model aimed at the description of instabilities due to load 
congestion that takes into account the topological proper- 
ties of SF networks. The results obtained point out that 
the network can freely handle traffic up to some critical 
average load {£)],. Above this level the network faces par- 
tial congestions that start to build up local bottlenecks in 
various places and small instabilities might trigger macro- 
scopic outages with a finite probability. Above a critical 
load value (i) 1 / ~ 0.82 any small instability leads to the 
whole network collapse. In the intermediate region of 
network load, the number of simultaneous line casualties 
follows a power-law resembling what has been observed 
in experimental studies of the Internet. We hope that 
our work will provide hints for accurate modeling of the 
Internet and the WWW large-scale traffic behavior. 
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